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Abstract 

Excited states are stationary localized solutions of the Gross-Pitaevskii equation with a 
harmonic potential and a repulsive nonlinear term that have zeros on a real axis. Existence and 
asymptotic properties of excited states are considered in the semi-classical (Thomas-Fermi) 
limit. Using the method of Lyapunov-Schmidt reductions and the known properties of the 
ground state in the Thomas-Fermi limit, we show that excited states can be approximated by 
a product of dark solitons (localized waves of the defocusing nonlinear Schrodinger equation 
with nonzero boundary conditions) and the ground state. The dark solitons are centered at 
the equilibrium points where a balance between the actions of the harmonic potential and the 
tail-to-tail interaction potential is achieved. 

1 Introduction 

The defocusing nonlinear Schrodinger equation is derived in the mean-field approximation to 
model Bose-Einstein condensates with repulsive inter-atomic interactions between atoms. This 
equation is referred in this context to as the Gross-Pitaevskii equation [9]. When the Bose- 
Einstein condensate is trapped by a magnetic field, the Gross-Pitaevskii equation has a harmonic 
potential. In the strongly nonlinear limit, referred to as the Thomas-Fermi limit pfl[TT], the Bose- 
Einstein condensate is a nearly compact cloud, which may contain localized dips of the atomic 
density. The nearly compact cloud is modeled by the ground state of the Gross-Pitaevskii 
equation, whereas the localized dips are modeled by the excited states. Asymptotic properties of 
the stationary excited states in the Thomas-Fermi limit are analyzed in this article. 

The Gross-Pitaevskii equation with a harmonic potential and a repulsive nonlinear term can 
be rewritten in the form 

ieu t + e 2 u xx + (1 — x 2 — \u\ 2 )u = 0, (1) 

where e > is a small parameter to model the Thomas-Fermi asymptotic regime. Let r] e be the 
real positive solution of the stationary equation 

e 2 r]"(x) + (1 - x 2 - rf £ {x))r] £ (x) =0, i£l. (2) 

Main results of Ignat & Millot [BJ [7] and Gallo Sz Pelinovsky 3] state that for any sufficiently 
small e > there exists a unique smooth positive solution rj £ £ C°°(R) that decays to zero as 
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\x\ — > oo faster than any exponential function. The ground state converges pointwise as e — ► 
to the compact Thomas-Fermi cloud 

/ \ ,. / \ / (1 - x 2 ) 1 / 2 , for \x\ < 1, 
%(x):=hmife(x) = | Q) for | x | >L ( 3 ) 

The ground state and the convergence of % to 770 is characterized by the following properties: 
PI < T) £ (x) ^ 1 for any i£i 

P2 For any small e > and any compact subset K C (—1,1), there is Ck > such that 

11% -Vo\\ci(K) < C K e 2 . (4) 
P3 For any small e > 0, there is C > such that 

\\Ve ~ VoUoo ^ Ce 1 / 3 , H\\ L ~ < Ce" 1 / 3 . (5) 

P4 There is C > such that ^ Ce 1 / 3 for any \x\ ^ 1 + e 2 ^ 3 . 

Properties [PI] and [P2] follow from Proposition 2.1 in [6]. Properties [P3] and [P4] follow 
from Theorem 1 in [3]. To clarify the proof of bound ([5]), we represent the ground state rj E {x) in 
the equivalent form 

r, £ (x) = e 1 ' 3 v e {y), V = ^r, (6) 

where v e {y) solves 

4(1 - e 2 / 3 y y:(y) - 2e 2 l 3 v' £ {y) + yv e {y) - v 3 e {y) = 0, y G (-ex,, e" 2 / 3 ). 
Let fo (y) be the unique solution of the Painleve-II equation 

H(y) + vMy) - "$(y) = o, yeR, 

such that 1/0 (y) = y 1 ^ 2 + C(y _1 ) as y —* 00 and fo(y) decays to zero as y — » — 00 faster than 
any exponential function. By Theorem 1 in [3J, v e is a C°° function on (— 00, e~ 2 / 3 ], which is 
expanded into the asymptotic series for any fixed N ^ 0: 

N 

My) = Y,z 2n/3 My) + e( 2N+1)/3 R N Ayl (?) 

n=0 

where {u n }^ =l are uniquely defined e-independent C°° functions on R and i?Ar ;£ (y) is the remain- 
der term on (—00, e~ 2 / 3 ]. It was proved in [3] that Un jE (z) = R^^ E {e~ 2 / 3 — e 2 / 3 z 2 ) is uniformly 
bounded for small e > in ff 2 (R)-norm. If we denote un, e (x) = Un, e (£~ 2 ^ 3 x) = RN, E (y), then 
the above arguments shows that there is Cat > such that 

H^V.elli 00 ^ Cn, \\ u N e\\L°° ^ Cjv£ -2 / 3 . 
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For any fixed N ^ 0, it follows from the above bounds that the remainder term e^ 2N+1 ^ 3 UN, £ ( x ) 
is smaller in C 1 (R) norm than the leading-order term uq(x) = ^o( e ~ 2 ^ 3 ~~ e~ 2 / 3 x 2 ). The error 
estimate © follows from ©, ©, and the fact that sup yeK + |fo(y) — V l ^ 2 \ < °o. 

We shall consider excited states of the Gross-Pitaevskii equation ([I]), which are real non- 
positive solutions of the stationary equation 

e 2 u"(x) + (l-x 2 -u 2 (x))u e (x) = 0, xeR. (8) 

We classify the excited states by the number m of zeros of u e (x) on M. A unique solution with m 
zeros exists near e = e m for e < e m by the local bifurcation theory [8], where e m is computed from 
the linear theory as e m = 1+ 1 2m , m £ N. Because of the symmetry of the harmonic potential, the 
m-th excited state is even on M. for even m£fJ and odd on R for odd m £ N. 

This paper continues the previous research on the ground state in the Thomas-Fermi limit 
that was developed by Gallo & Pelinovsky in [21 [3]. We focus now on the existence and asymptotic 
properties of the excited states as e —* 0. Using the method of Lyapunov-Schmidt reductions, 
we show that the m-th excited state is approximated by a product of m dark solitons (localized 
waves of the defocusing nonlinear Schrodinger equation with nonzero boundary conditions) and 
the ground state r/ E . The dark solitons are centered at the equilibrium points where a balance 
between the actions of the harmonic potential and the tail-to-tail interaction potential is achieved. 

Note that this paper gives a rigorous justification of the variational approximations found by 
Coles et al. in [1], where the m-th excited states was approximated by a variational ansatz in 
the form of a product of m dark solitons with time-dependent parameters and the ground state. 
Time- evolution equations for the parameters of the variational ansatz were found from the Euler- 
Lagrange equations. Critical points of these equations give approximations of the equilibrium 
positions of the dark solitons relative to the center of the harmonic potential and to each others, 
whereas the linearization around the critical points give the frequencies of oscillations of dark 
solitons near such equilibrium positions. Variational approximations were found in pQ to be in 
excellent agreement with numerical solutions of the stationary equation (jHj). 

This article is organized as follows. The first excited state centered at x = is considered 
in Section 2. Although existence of this solution can be established from the calculus of varia- 
tions, we develop the fixed-point iteration scheme to study this solution as e — > 0. The second 
excited state is approximated in Section 3. We will work with the method of Lyapunov-Schmidt 
reductions to find the equilibrium position of two dark solitons as e — > 0. Section 4 discusses the 
existence results for the general m-th excited state with m ^ 2. 

Before we proceed with main results, let us discuss some notations. If A and B are two 
quantities depending on a parameter e in a set £ , the notation A(e) = 0{B{e)) as e — > indicates 
that A(e)/B(e) remains bounded as e — ► 0. If A(x,e) depends oni£l and s E £ , the notation 
A(-,e) = Ol°°{B(e)) as e — > indicates that \\A{-,e)\\l°° / B{e) remains bounded as e — » 0. 
Different constants are denoted with the same symbol C if they can be chosen independently of 
the small parameter e. 
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2 First excited state 



The first excited state is an odd solution of the stationary equation ([8]) such that 

u £ (0) = 0, uJx) > for all x > 0, and lim uJx) = 0. (9) 

x— >oo 

Variational theory can be used to prove existence of this solution, similar to the analysis of Ignat 
&: Millot in [7]. Since we are interested in asymptotic properties of the first excited state as e — * 0, 
we will obtain both existence and convergence results from the fixed-point arguments. Our main 
result is the following theorem. 

Theorem 1 For sufficiently small e > 0, there exists a unique solution u e G C°°(M) with prop- 
erties (G|) and there is C > such that 

< Ce 2 / 3 . (10) 

L°° 

In particular, the solution converges pointwise as e — * to 

uo(x) := limu e (x) = rjo(x)siga(x) , 

e^0 



u £ — r/ e tanh 




Remark 1 Function v e (x) = tanh (^~^^J ^ s termed as the dark soliton. It is a solution of the 
second-order equation 

e 2 v"(x) + (1 - v 2 £ {x))v £ {x) = 0, x G R, 
which arises in the context of the defocusing nonlinear Schrddinger equation. 

The proof of Theorem Q] consists of six steps. 

Step 1: Decomposition. Let us substitute u e (x) = r] £ (x) tanh (^^^j + w e( x ) to the sta- 
tionary equation ([8]) and obtain an equivalent problem for w £ written in the operator form 

L £ w £ = H £ + N £ (w £ ), (11) 

where 

2a2 ,2 r , ,J/„u„„l,2 



L e := -e l d l x + x - 1 + 3r/ 2 (x) tanh 



.V2sJ ' 

H £ (x) := %(x) (r] 2 (x) — l) sech 2 ^— jL-^ tanh ^— ^=-^ + v / 2e?/e(2;)sech 2 



V2e 



and 



N e (w e )(x) = — 3r/ £ (x)tanh ( — j=- ) w £ (x) —w £ (x). 



Let x = y/2ez, where z G R is a new variable, and denote 

r? e (z) := r? £ (\/2ez), t& e («) := w £ (V2ez), H £ (z) := H £ (V2ez), N £ (w £ )(z) := N £ (w £ )(V2ez). 
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Step 2: Linear estimates. In new variables, operator L £ becomes 

L £ = --d 2 z + 2e 2 z 2 - 1 + 3rj 2 (z) tanh 2 (z) = L + U £ (z), 

where 

L ■= _I^ + 2-3sech 2 (z) 

and 

U £ (z) := 2e 2 z 2 + 3(fj 2 (z) - 1) tanh 2 (z). 

Operator Lq is well known in the linearization of the defocusing NLS equation at the dark 
soliton. The spectrum of Lq in L 2 (M.) consists of two eigenvalues at and | with eigenfunctions 
sech 2 (z) and tanh(z)sech(z) and the continuous spectrum on [2,oo). For any / G L 2 dd (]R), there 
exists a unique Lq 1 f G ff 2 dd (IR) such that 

3O0: V/GL 2 dd (M) : \\L^f\\ H , < C\\f\\». (12) 

Let us consider functions that decay to zero as |z| —* oo with a fixed exponential decay rate 
a > 0. Let L^(R) be the exponentially weighted space with the supremum norm 

[|«>e||l£> := [|e a l'l«) E ||i<x>. 

The unique solution Lq 1 / for any / G L 2 dd (]R) is expressed explicitly by the integral formula 

Lq X }{z) = -2sech 2 (z) J\osh\z') ( J" /(z ,, )sech 2 (z ,/ )dz ,/ ^ dz' . 

For any fixed a > 0, it follows from the integral representation that the solution Lq 1 / decays 
exponentially with the same rate as / so that 

3C>0: V/GL 2 dd (M)nL-(IR): H^VlU? < C\\f\\ L ~. (13) 

Figure 1 shows the confining potential V e (x) = x 2 — 1 + 2>r] 2 {x) tanh 2 (z) of operator L e = 
—£ 2 d 2 + V £ (x) (solid line) and the bounded potential Vq(x) = — 1 + 3tanh 2 (z) of operator L = 
—e 2 d 2 + Vq(x) (dots) versus x. The confining potential V £ {x) has two wells near x = ±1 and a 
deeper central well near x = 0. The two wells near x = ±1 are absent in the potential Vq(x). 

Because of the confining potential, the spectrum of L £ is purely discrete (Theorem 10.7 in 
[5j). It contains small eigenvalues that correspond to eigenfunctions localized in the central well 
near z = and in the two smaller wells near z = 

We note that a similar operator at the ground state e £ 

L £ = -e 2 d 2 + x 2 - 1 + 3r/ 2 (x) 

was studied by Gallo & Pelinovsky [3], where it was shown that V £ (x) = x 2 — 1 + 3r/ 2 (x) > for 
all By property (P4), V £ (x) is bounded away from zero near x = ±1 by the constant of 

the order of 0(e 2 / 3 ). As a consequence, the purely discrete spectrum of L £ in L 2 dd (IR) includes 
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Figure 1: Potentials of operators L £ (solid line) and Lq (dots) for the first excited state. 

small positive eigenvalues of the order 0(e 2 / 3 ) with the eigenfunctions localized in the two wells 
near x = ±1 (see Theorem 2 in [3]). 

Thanks to the proximity of tanh 2 (z) to 1 near z = with an exponential accuracy in e, 

the potential V e {x) is similar to V e {x) near x = ±1 and satisfies for any fixed xq > 0: 

3C > : V(x) > Ce 2/3 , |x| > x . 
On the other hand, for any fixed zq > 0, property (P2) implies that 

3C > : sup \U £ (z)\ ^ Ce 2 . 

\z\^z Q 

Thanks to the positivity of V e {x) near x = ±1 and the proximity of the central well near x = 
in the potentials 14 (x) and Vq(x), the quantum tunneling theory [5] implies that the simple zero 
eigenvalue of Lq persists as a small eigenvalue of L e . This eigenvalue of L £ corresponds to an 
even eigenfunction. The other eigenvalue of Lq corresponding to an odd eigenfunction is bounded 
away from zero. 

All other eigenvalues of L e are small positive of the size 0(e 2//3 ). As a result, operator L e is 
still invertible on £ 2 dd (M) but bound (p2D is now replaced by 

3C>0: V/GL 2 dd (M): ||4~7lliP < C £ - 2 / 3 ||/|| i2 . (14) 

Note that the function L^ 1 f £ H^ dd (M.) has peaks near points z = i— ^; and z = 0. 

Step 3: Bounds on the inhomogeneous and nonlinear terms. By symmetries, we note 
that 

H £ eL 2 odd (R) and N £ (w £ ) : H 2 dd (R) ^ L 2 odd (R). 
We will show that for small e > and fixed a G (0, 2) there is C > such that 

||4lbni- < Ce 2/3 - (15) 
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Using the triangle inequality, we obtain 

ll^elb < ll%IU-ll(l-€)sech 2 (-)|| L 2 + ^e||^|| L ^||sech 2 (-)|| L 2. 
By properties (PI) and (P2), for small e > and fixed a G (0,2) the first term is estimated by 
||(1 - r) 2 )sech 2 (.)|| L2 < ||(1 - r? 2 )sech 2 (.)|| i2( |,|^-i/3 ) + ||(1 - r/ 2 )sech 2 (.)|| i2( | 2 |^-i/3 ) 

< ||1 - ^ 2 || L oo ( | x | <v ^ £2 /3 ) ||sech 2 (-)|| L 2 + a- 1 / 2 e- Q£ " 1/3 ||sech 2 (-)||L- 

By property (P3), the second term is estimated by Ce 2 / 3 . As a result, for any small e > there 
is C > such that ||^ £ || L 2 ^ Ce 2 / 3 . By similar arguments, H £ G L™(R) for any a G (0,2) and 
there is C > such that H-HUIl^ < Ce 2/3 . 

II c 1 1 ^ a 

To deal with the nonlinear terms, we recall that H 2 (M.) is Banach algebra with respect to 
multiplication in the sense that 

Vu,t> G H 2 (R) : \\uv\\jj2 ^ ||m||#2||u||#2 

For any w £ G H 2 (M.), we have 

\We{w £ )\\ L 2 < 3||?7 e ||Loo||t£; 2 || H 2 + \\w*\\ H 2 < 3||u) E ||^a + ||u) e |||-2. (16) 

Similarly, L£?(R) is a Banach algebra with respect to multiplication for any a ^ 0. 

Step 4: Normal-form transformations. Because we are going to lose e 2//3 as a result 
of bound (|14|) . we need to perform transformations of solution w e , usually referred to as the 
normal- form transformations. We need two normal- form transformations to ensure that the 
resulting operator of a fixed-point equation is a contraction. 

Let 

w £ = wi + w 2 + (f>ei w\ = Lq H e , w 2 = — 3Lq 1 f) £ tanh(z)ti) 2 . 
The remainder term (p e solves the new problem 

c E E = n £ +K(0e), (17) 

where the new linear operator is 

C e := L E + AU £ (z), AU e (z) := 6r/ e tanh(z)(i&i + w 2 ) + 3(i&i + w 2 f , 
the new source term is 

H £ := -U £ (wi + w 2 ) - 3% t&nh(z)(2wiW2 + w 2 ) - (wi + w 2 ) 3 , 
and the new nonlinear function is 

M e (<p £ ) ■= -3fi £ t&nh(z)(p 2 - 3(wi + w 2 )ip 2 - cp 3 .. 
Thanks to bounds (USD, and CE}, we have w 1} w 2 G # 2 dd (M) n L~(M) for fixed a G (0, 2) 

and 

3C>0: ||u>i||tf2 n£ «> < Ce 2/3 , \\w 2 \\ H 2 nL oo < (7e 4/3 . (18) 
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As a result, for any small e > 0, there is C > such that 

\\fi £ t&n\i(z)(2wiW2 + wl)\\ L 2 < Ce 2 , [|(t& a +w 2 ) 3 ||l2 < Ce 2 . 

Let us now estimate the term U £ (w\ + W2) in £ 2 (K). By properties (PI) and (P2), for small e > 
and fixed a £ (0, 2) there are constants C(a), C(a) > such that 



1 1 U £ Wj\\ L 2 < 2e 2 ||z 2 t() J -|| i 2 + 3||(i7 2 - l)wj\\ L 2 

< e 2 C(a) 11^- || L go +3||(1 -€H-|Il2(|*|^ £ -V3) + 3||(1 - €K'IIl2(M> £ -i/3) 

^ e 2 C(a) 11^- || L oo +3||1 -^ILoo ( | a .| < ^ e a/a ) ||u)j||La + 3or 1 / 2 e- ae ~ 1/ *\\w j \\ L oo 

< C , (a)e 4 / 3 ||T& i || i a nJ;g ,, j = 1,2. 

In view of bound (JTHJ) , for any small e > there is C > such that 

\\U £ (wi + w 2 )\\ L 2 ^Ce 2 . (19) 

Combining all together, we have established that TL e € L 2 dd (M) and for any small e > 0, there is 
C > such that 

||H e || L2 ^ Ce 2 . (20) 

For the nonlinear term, we still have N £ {<p £ ) : H 2 dd (R) ^ L 2 dd (M). Thanks to bound (|T8"|) . 
for any <^ e G Bs(H 2 dd ) in the ball of radius 5 > 0, for any small e > there is C(5) > such that 

\\N £ (0 £ )\\ L ^C(8)\\0 £ f H2 . (21) 

Similarly, we obtain that A4 is Lipschitz continuous in the ball Bs(H 2 dd ) and for any small e > 
there is C(<5) > such that 

V<£ £ , </? e G B$(H 2 dd ) : ||A4(^)-A4(^)|| L 2 <C(<5) (11^11^2 + 11^11^)11^-^11^2. (22) 



Step 5: Fixed-point arguments. Thanks to bound (|18p and Sobolev embedding of if 2 
to L°°(R), |AC/ e (z)| is as small as 0(e 2 / 3 ) in the central well near z = and is exponentially 
small in e in the two wells near z = As a result, small positive eigenvalues of L £ of the 

size 0(e 2//3 ) persist in the spectrum of C £ and have the same size, so that bound (|14|) extends to 
operator L £ in the form 

BOO: V/GL 2 dd (M): \\CT £ x j\\n^ C^ 2/3 ||/ll^- (23) 

Let us rewrite equation (fTTf) as the fixed-point problem 

£ e F 2 dd (M) : & = + £-W £ fe). (24) 

The map yj £ 1— » C~ 1 J\f e (<p e ) is Lipschitz continuous in the neighborhood of G H 2 dd {U). Thanks 
to bounds (f2Tj) and ([23]) . the map is a contraction in the ball Bs(H 2 dd ) if 5 -C e 2 / 3 . On the 
other hand, thanks to bounds (|20p and (j23h . the source term C~ l TL £ is as small as C(e 4 / 3 ) m 
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L 2 norm. By Banach's Fixed-Point Theorem in the ball B$(H 2 dd ) with 5 ~ e 4//3 , there exists a 
unique ip e G if 2 dd (R) of the fixed-point problem (I24p such that 

3C>0: \\<p e \\H2 < Ce 4/3 . 
By Sobolev's embedding of if 2 (R) to C 1 (M), for any small e > there is C > such that 

||w £ ||l°° = ||^ £ ||i°° < C\\wi + w 2 + ^ £ ||h 2 < Ce 2 / 3 , 

which completes the proof of bound (|10p . 

Step 6: Properties ([9]). Solution w e constructed in Step (5) is a odd continuously differ- 
entiable function of z on R vanishing at infinity, so that u e (0) = and lim x ^ 00 u e (x) = 0. By 
bootstrapping arguments for the stationary equation ([5]), we have u £ G C°°(R). It remains to 
prove that u £ (x) is positive for all x G R + . 

Recall that r/ e (x) > for all By property (P4) and bound (fTUj) . there is C > such 

that u e (x) ^ Ce 1/3 for all x G [1, \J 1 + e 2 / 3 ]. We shall prove that u E (x) > for all x > \A + e 2 / 3 . 
Assume by contradiction that there is xo > yl + e^ such that u e (xo) = and u' e (xo) < 0. (If 
u' £ (xo) = 0, then u e (x) = is the only solution of the second-order equation (JSJ.) The continuity 
of u £ (x) implies that u e (x) < for every x G (xo,xo) for some xo > ^o- Using the differential 
equation {SJ, we obtain 

u' £ (x) = -^(x 2 - 1 + u 2 £ {x))u £ (x) < 0, x G (x ,xo). 

Then, u' £ (x) ^ ^(xo) < 0, so that u £ (x) is a negative, decreasing function of x for all x > xo 
with xo = oo. This fact is a contradiction with the decay of u £ (x) to zero as x — ► oo. Therefore, 
u £ {x) > for all x G E + . 

Combining results in Steps (5) and (6), we conclude that u £ (x) is the first excited state of the 
stationary equation (|8|) that satisfies properties Q. 



3 Second excited state 

The second excited state is an even solution of the stationary equation ([8]) such that 

u £ (x) > for all |x| > xo, uJx) < for all |x| < xq, and lim u £ (x) = 0. (25) 

x— >oo 

Here xo > determines a location of two symmetric zeros of u £ (x) at x = ±xo The second excited 
state is approximated as e — ► by a product of two copies of dark solitons (Remark [[]) placed 
at x = ±a with a ~ xo as e — > 0. Our analysis is based on the method of Lyapunov-Schmidt 
reductions, which gives existence and convergence properties for the second excited state, as well 
as an analytical expansion of a for small e. 



Theorem 2 For sufficiently small e > 0, there exists a unique solution u £ G C°°(R) with prop- 
erties A25\) and there exist a > and C > such that 



rj E tanh 



V2e 



tanh 



• + q 



< Ce 2 / 3 



(26) 
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and 

e / , , 1 , 3 



o = --^(log(e) + -log|log(e)|--log(2)+o(l)J as e - 0. (27) 



Remark 2 Since a — ► as e — > w;Me r/ e (x) ~ 1 near x = 0, we /iaue 

x = a + C(e 5/3 ) as e 0. 



Remark 3 Exactly the same asymptotic expansion Jl?7| ) has been obtained with the use of the 
averaged Lagrangian approximation and has been confirmed numerically fl^. 

The proof of Theorem [2] follows the same steps as the proof of Theorem [1] with an additional 
step on the Lyapunov-Schmidt bifurcation equation. 

Step 1: Decomposition. Let a £ (0,1) and substitute 

u £ (x) = ij £ (x) tanh f ^_ ~\ tanh f ~\ + w £ (x) 

to the stationary equation ([8]). The equivalent problem for w £ takes the operator form 

L £ w £ = H £ + N £ (w £ ), (28) 



where 



L £ := — e d x + x — 1 + 3i] £ (x) tanh (z+) tanh (z~), 

H £ := rj e (x)(r] e (x) — 1) tanh(z + ) tanh(z_) (sech 2 (z + ) + sech 2 (z_)) 
+n e (x)sech 2 (2;_|_)sech 2 (2;_) (l — r/ 2 (x) tanh(z + ) tanh(z_)) 
+\f2er] £ {x) (tanh(z + )sech 2 (z_) + tanh(z_)sech 2 (z + )) , 



and 



N £ (w £ ) = — 3rf e (x) tanh(z + ) t&nh(z-)w £ (x) — w £ (x) 
with the following notations 



x 



z± = z±(, z = — Er , C 



y/2e' V2e' 

We again denote the functions in z by hats. We shall assume a priori that 

3/9 €(0,1): a^V2/?e 2 / 3 , 
3C>0: e" 2 ^ Ce|log(e)| 1 / 2 . 

Note that bounds (|29|) imply that a — ► and ( — ► oo as e — > 0. 



(29) 
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Step 2: Linear estimates. In new variables, operator L £ becomes 

L e = -^d 2 + 2e 2 z 2 - 1 + 3fj 2 £ {z) tanh 2 (z + Q tanh 2 (z - C) = L (C) + U e (z, C), 

where 

Lo(C) = + 2 - 3sech 2 (z + C) - 3sech 2 (z - C) 

and 

U e (z, C) = 2e 2 z 2 + 3(fi 2 (z) - 1) tanh 2 (z + () tanh 2 (z - C) + 3sech 2 (z + C)sech 2 (z - (). 

Operator £o(C) has now two eigenvalues in the neighborhood of for large £ because of the 
double-well potential centered at z = ±£. If C is large, the geometric splitting theory [10] implies 
that the eigenfunctions tp (z) of operator Lo(C) corresponding to the two smallest eigenvalues 
are given asymptotically by 

?±( \ ^(z - () ± tpo(z + () ,„ , _2C\ a fon\ 

Wq {z) = ■= hC L oo(e ^) as ( ^ oo, (30) 

v2 

where V>o(-2) = -^sech 2 (2i) is the L 2 -normalized eigenfunction of Lq = —\d\ + 2 — 3sech 2 (z) for 
the zero eigenvalue. 

Note that ^q(z) is even and ^q(z) is odd on R. For the second excited state, we are looking 
for an even solution w £ {z). Since a is not specified yet, we add the condition (^>jj" ,w e ) = and 
define a constrained subspace of -ff 2 ven (R) by 

X = {w £ G H 2 vcn (R) : w e ) = 0}. 

Let Po be an orthogonal projection operator to the complement of in L 2 vcn (R). Since eigen- 
function x/jq is odd and the rest of spectrum of -Lo(C) is bounded from zero, for any / G L 2 ven (R), 
there exists a unique PqLq 1 (()Pof G H 2 ven (M.) such that 

3C> : V/ G L 2 vcn (R) : \\P L^ HO^o/ll^ < C||/|| L2 . (31) 

Let us consider functions that decay to zero as \z — Q,\z + Q — > oo with a fixed exponential 
decay rate a > 0. Let L^(R) be the exponentially weighted space with the supremum norm 

\\w E \\ L oo ■= sup e ailz -^\w £ (z)\ + sup e a(|z+cl) |w £ (z)|. 

For fixed a > and £ > 0, the unique solution PqLq 1 (()Pof decays exponentially with the same 
rate as / so that 

3C> : V/ G Ll cn (R) n L~(M) : ||P ^ HCWlli^ < C|I/IU- C - (32) 

Figure 2 shows the potential V £ (x) = x 2 — 1 + 3r/ 2 (x) tanh 2 (z + C)tanh 2 (z — £) of operator 
L e = -e 2 <9 2 + K(x) (solid line) and the potential V (x) = 2- 3sech 2 (z + £) — 3sech 2 (z - Q of 
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Figure 2: Potential of operator L e (solid line) and Lq (dots) for the second excited state. 

operator Lq = —£ 2 d 2 + Vq[x) (dots) versus x. The bounded potential Vq{x) has two wells near 
x = ±a, whereas the confining potential V e (x) has four wells near x = ±a and x = ±1. 

Again, the spectrum of operator L e with a confining potential is purely discrete. The two 
wells of the confining potential V £ (x) near x = ±1 are 0(e 2 / 3 )-close to zero but still positive 
thanks to property (P4) and the fact that tanh(z ± £) = 1 with exponential accuracy in e if 
£ ^ j3e~ 1 /^ for fixed (3 € (0, 1). Therefore, for any fixed xq > 0, we have 

3C > : V(x) ^ Ce 2/3 , \x\ ^ x . (33) 

On the other hand, by property (P2) for any £ £ (0,e™ 2 / 3 ), we have 

3C > : sup \U £ (z, C)| C(e 2/3 + e" 4C ). (34) 

Thanks to properties (i33|) and (i34"|) . the quantum tunneling theory [5] implies that the two small 
eigenvalues of Lq persist as two small eigenvalues of L £ with two eigenfunctions ip^ that satisfy 
asymptotically 

$f(z) = 4>q{z) +0 L oo(e 2 / 3 ) as e -> 0, (35) 

thanks to a priori bound (|29p and the exponential smallness of tp^iz) in e near z = 

Let P e be an orthogonal projection operator to the complement of ip^~ in L 2 ven (R). Because 
of the small 0(e 2 ^ 3 ) eigenvalues of L e , bound (|31h is now replaced by 

30 0: V/aL(K): \\P £ L- l P £ f\\ H ^ Ce~ 2 l 3 \\f\\ L ,. (36) 

The function P £ L~ 1 P £ f G H 2 ven (M>) has peaks in all four wells near points z = and 2 = ±C- 

Step 3: Bounds on the inhomogeneous and nonlinear terms. From the symmetry of 
terms in H £ and N e (w £ ), we have 

N £ (w £ ) : H 2 vcn (R) ' * L 2 cvcn (R) and JT £ G ^ c 2 vcn (M). 
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Under a priori bound (|29p . we first show that there is C > such that 

\\H £ \\ L 2 ^ Ce 2 ' 3 . (37) 

The upper bound for the first term in H e involves estimates of 

h{z) := (1 - 77 e 2 (z))(sech 2 (z + C) + sech 2 (z - C)), 

which may create a problem since £ — > oo as e — > and ?) £ (z) — > as |z| — >• oo. By properties 
(PI) and (P2), for any a G (0,2), £ < /3e -1 / 3 for any G (0, 1), and any small e > 0, there is 
constant C > such that 

II-^iIIl 2 ^ \\h IIl2(| z |< £ -i/3) + \\h 11x2(1^-1/3) 

< II 1 -^llLoo ( | x .| <v ^ e 2/ 3) ||sech 2 (2; + ) +sech 2 (z_)|| L 2 

+ a-Vh- a l e ~ 1/a -Q\\sech 2 (z + ) + sech 2 (z_)|| L - c 

^ Ce^ 3 . 

Thus, the condition £ ^ (ie~ 1 ^ 3 from a priori bound (|29p is sufficient to keep I\ small in L 2 . 
The upper bound for the second term in H e involves the estimate of the overlapping term 

h(z) '■= sech 2 (z_|_)sech 2 (2:_). 

Under a priori bound (|29p . this term is estimated by 

1/2 



\h\\L 2 ^ 



Jr J 

/oo \ 1/2 

sech 4 (u)sech 4 (u + 2Q)duj < Ce~ 2<: ^ Ce\ Log(e> 



The last term in H £ is proportional to er}' £ and is handled with property (P3) to give (|37p . By 
similar arguments, H e G L^(M) for any a G (0, 2) and £ ^ f3e~ 1 ^ 3 for any [5 G (0, 1) and for any 
small e > there is C > such that ||.H e ||.L«> ^ Ce 2 / 3 . 

The nonlinear terms in N e (w e ) are handled with the Banach algebra of H 2 (M.), so we obtain 

\\N E (w £ )\\ L 2 ^311^11^ + 11^11^. (38) 

Step 4: Normal-form transformations. Unlike step (4) in the proof of Theorem (H we 
need to perform a sequence of two normal-form transformations because the orthogonal projection 
operator to the one-dimensional subspace spanned by an even eigenfunction for the smallest 
eigenvalue of Lq has to be changed to the projection operator associated with an eigenfunction 
of a new linearization operator. For the sake of short notations, we combine both normal-form 
transformations and write them together. 

Let w e = w\ + u>2 + ¥>e with 
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where 

G £ := -3f/ £ tanh(z + £) tanh(z - C)u>? + (V e ~ Pq)H £ 

and V e is a new orthogonal projection operator introduced below. 
The remainder term e solves the new problem 

C £ E = H E +K(0e)+S E , (39) 

where the new linear operator is 

C £ := L £ + AU £ (z), AU £ (z) := 6%tanh(z + Q tanh(z - ()(wi + w 2 ) +3(u>i + w 2 ) 2 , 
the new source term is 

H £ := -U £ {w l +w 2 ) - 3%tanh(z + C)tanh(z -C)(2wiw 2 + - (u>i + w 2 f + (V £ - P )G £ , 

the new nonlinear function is 

M £ (ip £ ) := -3% tanh(z + Q tanh(z - ()0 2 - 3(u>i + w 2 )0 2 - 0^, 

and the new one-dimensional projection is 

S £ :={I-V £ ){H £ +G £ ). 

If t&i, u> 2 € H 2 (R) n L™ C (R) satisfy bounds <@2J) below, then AU £ (z) is as small as C(e 2 / 3 ) in 
the two wells near z = ±£ and is exponentially small in e in the two wells near z = i^g- Let ^jr 
be the eigenfunctions of C £ for the two eigenvalues continued from the two smallest eigenvalues 
of Lq. The proximity of the potential wells and expansion (|35|) imply that 

$f{z) = 4 ± (z) + LO c( e 2 / 3 ) = i>±(z) + Lac (e 2 / 3 ) as e -> 0. (40) 

Let be an orthogonal projection operator to the complement of ^+ in L 2 ven (K). Thanks to 
expansion (|4"0"]) . we have 

3C>0: ||P £ -P ||l2^ l2 ^Ce 2 / 3 . (41) 

Thanks to bounds ^3T]>, ([32]), (JHTJ), and flU}, we have w x ,w 2 G #evcnW n for any 

a G (0, 2) and C «S /3e~ 1/3 for any /3 G (0, 1) such that 

3C>0: ll^iH^nico <Ce 2 / 3 , ||to 2 || H2ni -^ Ce 4/3 . (42) 

As a result, for any small e > 0, there is C > such that 

11% tanh(z + C) tanh(z — ()(2wiw 2 + w 2 )\\ L 2 < Ce 2 , 

\\(wi + w 2 ) 3 \\ L 2 ^ Ce 2 , 
\\(V £ - P )G £ \\ L 2 < Ce 2 . 
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Let us now estimate the term U £ Wj in L 2 (M) for any j = {1,2}. By properties (PI) and (P2), 
for any a £ (0, 2) and £ ^ /3e -1 / 3 for any j3 £ (0, 1), and for any small e > 0, we have 



\e 2 z 2 w 



\ L 2 < Ce 2 C 2 ||^-|U- < C^/Sii^n 



Q,C " ■> " a, C ' 

a2\- II i ii/i -2\„~ 



||(1 - fl £ )wj\\ L 2 ^ ||(1 - %)«>il| £ a (w<e -i/S) + ||(1 - %)«y|lra(|*|>s-V3) 

^ 111 -^llL-(|x|<V2e2/3)ll^llL 2 +a- 1/2 e- Q ( £ " 1/3 -^||u; i || L oo c 

«S Ce 4/3 ||^llL2nL- c , 

||sech 2 (z_ ) _)sech 2 (z_)t£) £ || i 2 ^ Ce _4 ^||ti) £ || i 2 ^ Ce 2 \ log(e)|||w e || £ 2. 

In view of bound (|42p . for any small e > 0, there is C > such that 

||t7 £ (^i+u> 2 )|| L 2 ^Ce 2 . (43) 

Combining all together, we have established that Tt £ £ L 2 ven (M) and for any small e > 0, there 
is C > such that 

\\H £ \\ L 2 < Ce 2 . (44) 

For the nonlinear term, we still have Af e (<p e ) : H 2 ven (R) i-> L 2 ven (M). Thanks to bound (02]), 
for any <^ £ £ B$(H 2 vcn ) in the ball of radius 5 > and for any small e > 0, there is C(5) > such 
that 

\\K(0e)\\iP<C(8)\\&\\h- (45) 

Step 5: Fixed-point arguments. Because At/ e (z) is exponentially small in g near z = 
i^-j small positive eigenvalues of L £ of the size 0(e 2 / 3 ) persist in the spectrum of C £ and have 
the same size. As a result, bound (|36p extends to the operator C £ in the form 

3C> : V/ € LLnW : HTU^/Hfl* < Cg- 2 / 3 ||/|| x2 , (46) 



where the new projection operator V £ is used. As a result, we rewrite equation (139j) as the 
fixed-point problem 

£ £ #Ln W : ^ £ = VeC- x V e (H £ + A4(&) + S £ ) . (47) 
subject to the Lyapunov-Schmidt bifurcation equation 

T £ := (4 + ,(W £ +A4(^)+5 £ )) L 2 = $i,{H e + G e + H e +N e {<f> e )))v =0. (48) 



The map (p £ i— > V £ Cj 1 V £ M £ ((p £ ) is Lipschitz continuous in the neighborhood of £ H 2 ven , 
Thanks to bounds (|45p and (JM|), the map is a contraction in the ball Bs(H 2 ven ) if 5 <C e 2//3 . 
On the other hand, thanks to bounds (|44p and (|46p . the source term V £ C- X V £ H £ is as small as 
0(£ 4//3 ) in L 2 norm. Furthermore, V £ S £ = 0. 

By Banach's Fixed- Point Theorem in the ball B$(H 2 vcn ) with 5 ~ e 4 ' 3 , for any (a, £) satisfying 
a priori bounds (f29j) and sufficiently small e > 0, there exists a unique <y3 £ £ ^T| ven (R) of the fixed- 
point problem (j47j) and C > such that 



In* < Ce 4 / 3 . 



15 



By Sobolev's embedding of if 2 (R) to C 1 (M), for any small e > there is C > such that 

||w £ ||l°° = H^eIIl 00 s£ C\\wi + w 2 + fie\\H 2 < Ce 2/3 , 

which completes the proof of bound (|26p for any (a, satisfying a priori bounds (|29p . It remains 
to show that bounds (J29]) are satisfied by solutions of the Lyapunov-Schmidt bifurcation equation 



Step 6: Lyapunov Schmidt bifurcation equation. To consider solutions of the Lyapunov- 
Schmidt reduction equation, we rewrite ()48|) in the form 

where 

?P = (4+,(g £ + w £ +a4(^))) L 2 + (4 + -V'o + 1 ^)l- 

We will show that there exists a simple root of Te in a > 0, which satisfies the asymptotic 
expansion (12T|) and that this root persists with respect to the perturbations in T e . If a satisfies 



the asymptotic expansion ([27|) . then a = 0(e\ log(e)|) and e~ 2C - = 0(e\ log(e)\ 1 ^ 2 ) so that a priori 
bounds (f29|) are satisfied. 

For convenience, we recall (130)) and write 

K e = ^-^P = ((sech 2 (z+) +sech 2 (z_) + LOO (e~ 2C ),4) i2 . (49) 

In what follows, we compute the leading order of 1Z and account the error of the size OL°°(e~ 2< *) 
in the end of computations. From explicit definition of H £ , the leading-order part of lZ e is written 
in the form 

KP = ((sech 2 (z + ) + sech 2 (z_),i? £ ) L 2 

= / „«(*)(„?(*) - 1) tanh( 2+ ) ta„h( 2 _) (sech>( 2+ ) + sech V)) 2 

Jm 

+ / ?? £ (x)sech 2 (2;_|_)sech 2 (z_) (l — rj 2 (x) tanh(z + ) tanh(z_)) (sech 2 (z + ) + sech 2 (z_)) dz. 
Jm. 

After the change of variables u = z — £ = z_ = z + — 2£ and the use of symmetry on z 6 R, the 
first and second terms in 1Z £ give 



/1 + /2 



r] e (x)(r] 2 (x) - 1) tanh(u) tanh(u + 2£) (sech 2 (u) + sech 2 (u + 2()) du 

/oo 
rf e {x) (tanh(n)sech 2 (ii + 2Q) + tanh(u + 2C)sech 2 (u)) 



x (sech 2 («) + sech 2 (w + 2£)) du 



3\/2, 



;i + 7 ] 2 (x))7 1 ' £ (x)sech i (u) (l + L oo( e - 2 ^ du 
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where x = \^2e(u + £). Thanks to the exponential decay of sech 4 (w) and property (P3), we have 
3V2s f< 

i 

C 



h + h 



J (1 + 7? e 2 (x))r ? ;(x)sech 4 (n) (l + L oo (e~ 2 ^ du + 0(e 2 / 3 e~^). (50) 



2 

On the other hand, thanks to property (P2) for £ /fe -1 / 3 for any j3 G (0, 1), we have 

%(x) = l + L oc(e 4 / 3 ), t/ b (x) = -:r(l + L oc( e 4 / 3 )), Vs e [0, 2\/2eC]- 
As a result, we obtain 



/1 + /2 



-6e 2 y (C + -u)sech 4 (n) (l + Olo.^ 3 ,e~ 2C: )^ du + 0(e 2/3 e~^) 
-4V2ea ( 1 + 0(e^ 3 , e^)) + 0( e 2 /V 4 ^). 



Performing similar computations for the third term in lZ e gives 

I 3 = 2 J r? £ (x)sech 4 (u)sech 2 (n + 2C) (l -?? 2 (x)tanh(u)) (l + O (e~ 2C -)^ du 



-c 

32e" 4 ^ (l + 0(e A / 3 ,e-^))+0(e-^). 



Recalling now (|39j) . we have thus obtained that 

ft £ = -4V2ea ( 1 + 0{e 2 ' 3 , e^)) + 32 e - 2 ^ 1 ( 1 + 0(e A / 3 , e~ 2 ^ 



(2) 

Analyzing similarly the error coming from the other term Ti, in the Lyapunov-Schmidt reduction 



equation (j48J) , we rewrite this equation in the form 
2^/2 



^ £ = -AV2ea (l + 0(e 2 / 3 , e^)) + We- 2 ^' 1 (l + Q(e 2 / 3 , e~ 2 ^)) = 0. (51) 



Taking a natural logarithm of J- £ = 0, we obtain 

2^2a^ 1 + log(a) = - log(e) + \ log(2) + 0(e 2 / 3 , e"^ ). 
Let a = — ^elog(e)£/ and rewrite the problem for U: 

_ log(lQ _ log|log(e)| _ 31og(2) / 2/3 2C \ 

17 21og( £ )- 1+ 21og( £ ) 21og(.)l 1 + 0(£ ' 6 } J- (52) 



The remainder term is continuous with respect to e for small e > 0. There exists a root of (|52p 
at J7 = 1 for e = 0. By the Implicit Function Theorem applied to equation (i52"j) for small e > 0, 
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there exists a unique root U(e) such that U(e) is continuous in e > and lim £ ^ ^( e ) = 1- To 
estimate the remainder term, one can further decompose 

21og(e) 

and rewrite the problem for V: 

V log|log(e)| " log|log( £ )| l 1+0(£ ' e >)■ (53) 

Again, there is a root of (|53p at V = for e = 0. By the Implicit Function Theorem applied to 
equation (f53|) for small e > 0, there exists a unique root V(e) such that V(e) is continuous in 
e > and lim £ j^o V(e) = 0. As a result, for small e > there is a root of the nonlinear equation 
(|5ip . which admits the asymptotic expansion (|27p . 

Step 7: Properties (1251) . The uniform bound (|26p has again the order of 0(e 2 / 3 ). Using 
the same analysis as in Step 6 of the proof of Theorem [H we prove that u e (x) is strictly positive 
for any \x\ 1. Therefore, there exist only two zeros of u £ (x) on R and the two zeros x = ±xo 
are located near x = ±a (Remark [2]). Additionally, u £ G C 1 (M) and the bootstrapping arguments 
give u e G C°°(R). Combining all together, u £ (x) constructed above is the second excited state of 
the stationary equation ([8]) that satisfies property (|25p . 



4 Construction of the m-th excited state with m ^ 2 

The m-th excited state is constructed similarly to the proof of Theorem [2j The relevant decom- 
position is a product of m dark solitons and the ground state in the form 

m 



/ x — a \ 

u £ (x) = i] £ (x) Y\ tanh y ^_ 3 J + w £ (x), 



(54) 



i=i 

where parameters {aj}^L 1 are to be found from m constraints on the fixed-point problem for the 
remainder term w £ (x) . Assuming that all a,j are distinct and distributed according to the a priori 
bounds 

f 3(3 G (0, 1) : a, < ^/2> 2 / 3 , j G {1, 2, m} 
| 3C>0: e-^K+i-o^e 1 ^ C7e 2 |log(e)|, j G {1, 2, m — 1}, 

the relevant potential of the Schrodinger operator 

1 m ci 

L (ai, ...,0m) = --$1 + 2 - 3 V^sech 2 ^ - zj), Zj = —jL- 

has to wells and supports to eigenvalues in the neighborhood of 0. The m constraints follow from 
m projections to the corresponding eigenfunctions for the m smallest eigenvalues. Although the 
computations of these reductions are long and cumbersome, these computations are expected to 
recover the same leading order as the Euler-Lagrange equations obtained by Coles et al. PQ, 



18 



where only pairwise interactions contribute to the leading order. Asymptotic expansions of 
solutions of these equations are constructed in pQ and compared to the numerical approximations 
for m = 2 and m = 3. 

Spectral stability of the excited states in the limit e — > is also a physically important and 
mathematically interesting problem. Variational and numerical approximations in [I] suggest that 
the purely discrete spectrum of the spectral stability problem associated with the m-th excited 
state has a countable set of eigenvalues, which are close to eigenvalues associated with the ground 
state, and m additional pairs of eigenvalues. The m additional pairs are related to the Jacobian 
of the reductions equations (j55() : one pair remains bounded as e — > and (m — 1) pairs grow like 
log(e) as e — > 0. Unfortunately, the rigorous studies of the asymptotic properties of eigenvalues 
are difficult even for the linearization of the ground state [2j. Therefore, the characterization 
of asymptotic properties of eigenvalues associated with the excited states will remain an open 
problem for further studies. 
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